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. Abstract 

' A numerical investigation of a non-commutative field theory defined via the spectral 

action principle is conducted. The construction of this triple relies on an 8-dimensional 
Clifford algebra. Following to the standard procedure of non-commutative geometry, the 
spectral action is computed for the product of the triple {A^^'Hi^'Di) with a matrix- valued 
pi • spectral triple. Using Monte Carlo simulation we study various quantities such as the energy 

D , density, the specific heat density and some order parameters varying the matrix size and the 

'■^ ' independent parameters of the model. 



1 Introduction 

fSJ ' The main object of this work is a particular non-commutative field theory which is derived using 

the spectral action principle and then treated numerically. Non-commutativity can be found 
in many fields of physics like quantum field theories, string theory, condensed matter physics. 
The first application of non-commutativity into physics is dated from the middle of the last 
century inspired by the ideas of quantum mechanics, where starting from classical mechanics, the 
commutative algebra of functions on the phase space is replaced by a non-commutative operator 
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?H ' algebra on a Hilbert space. The duality between ordinary spaces M and proper commutative 

algebras is expressed by the Gel'fand-Naimark theorem which states the fact that the algebra of 
all continuous functions on M is the only possible type of commutative C*-algebra. Additionally, 
given a commutative C*-algebra C, it is possible to reconstruct a Hausdorff topological space M 
in order to obtain that C is the algebra of continuous functions on M. The study of commutative 
C*-algebras is equivalent to the study of topological Hausdorff spaces. The previous duality 
has inspired the identification, in non-commutative geometry, of some algebraical objects as a 
category of non-commutative topological spaces. Alain Connes [1], one of the founders of non- 
commutative geometry, has proposed a candidate for the objects of such category, the spectral 
triples {AyHjV) [2], composed by an algebra A, an Hilbert space Ti on which A is represented 
and an selfadjoint operator T>. In fact, every compact oriented Riemannian manifold can be 
used to define a spectral triple, this kind of manifold M characterizes a Dirac operator on 
self-adjoint Clifford module bundles over M. Connes, after a conjecture in 1996 f6] and some 
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considerable attempts of Rennie and Varilly proved the so called reconstruction theorem 
[3] for commutative spectral triples satisfying various axioms, showing that exists a compact 
oriented smooth manifold X such that A = C'^{X) is the algebra of smooth functions on 
X and every compact oriented smooth manifold emerges in this way. Pushed by the aim of 
reformulating the standard model of particles in a non-commutative way [4, 6j, Connes has 
introduced the almost-commutative spectral triple extending the axioms of the reconstruction 
theorem to a non-commutative algebra. A first attempt to formulate a field theory for a truly 
non- commutative algebra was obtained replacing in the usual field theory action the point- wise 
multiplication of the fields with a non-commutative one, namely a ^-product. The fields now 
belongs to Mq, a vector space defined by an enough regular class functions on Mg equipped with 
the Moyal product: 



Where is a skew-symmetric matrix. A very important question about non-commutative 
quantum field theory |S], is whether or not the quantum theory is well-defined or in other 
words if it is renormalizable or not. At first sight, the non- locality of the non-commutate 
action induced by the ★-product in the position space can induce us to fear some problems 
for the renormalization. In fact it was discovered [8] that after computing the Feynman rules 
for such theory and deriving the loop amplitude we find that the non-local interaction terms 
in the action induce an oscillatory factors (involving loop momenta) in the Feynman integrals. 
Studying the structure of Feynman diagrams for the action, Filk [13] has showed that are present 
two types of loop diagram the planar and non-planar diagrams. The planar diagrams do not 
have this oscillatory factors coming from the non-local interaction terms, and therefore the 
corresponding integrals are the same as in usual quantum field theory. On the other hand, all 
non-planar diagrams have the oscillatory factors involving loop momenta. Due to this terms 
the renormalization of quantum field theories on the non-commutative is not achieved and 
these models show a phenomenon called UV/IR- mixing [11] . Chepelev and Roiban [12] analyses 
UV/IR-mixing to all orders, the conclusion of the power-counting theorem is that field theories on 
non- commutative M" are not renormalizable if the divergence of their commutative counterparts 
are higher than logarithmic. 

A great step towards the non-commutative field theory was made when H.Grosse and 
R.Wulkenhaar [23], found a non-commutative 99^-theory renormalizable action which develops 
additional marginal coupling, corresponding to an harmonic oscillator potential for the real- 
valued free field if on Rq : 



Where x = 2Q~^ ■ x, A € M G [0,1], and ^ is a real parameter. Using the Moyal ma- 
trix base, which turns the ^-product into a standard (infinite) matrix product, H.Grosse and 
R.Wulkenhaar were able to prove the perturbative renormalizability of the theory [25]. After- 
ward, R.Wulkenhaar et al. [14] found an alternative simpler normalization proof using multi- 
scale analysis in matrix base, showing the equivalence of various renormalization schemes. A last, 
but useful, renormalization proof was formulated using Symanzik type hyperbolic polynomials 




f{x + ]^Q-k)g{x + y)e'^^'y^ 
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The non-commutative model treated in this work is a sort of extension, via spectral action 
principle, of the scalar W-G model, in which we are interested to formulate a Yang-Mills theory 
in renormalizable way on Moyal space. We can expect that usual Yang-Mills theory on Moyal 
space without modifications of the action by something similar to an oscillator potential, to be 
not renormalizable [11]. Additionally, the Moyal space with usual Dirac operator is a spectral 
triple, the corresponding spectral action was computed in [TTI, with the result that it is the 
usual not renormalizable action on Moyal plane. In [18] H.Grosse and R.Wulkenhaar, in order 
to obtain a gauge theory with an oscillator potential via the spectral action principle, used a 
Dirac operator constructed using the statement = H where the four dimensional Laplacian is 
substituted by the four dimensional oscillator Hamiltonian H = — A + ri^||x|p. The idea behind 
is that the spectral dimension is defined through the Dirac operator so the spectral dimension 
defined by such Dirac operator is related to the harmonic oscillator phase space dimension. 
It turns out that to write down an Dirac operator, so that its square equals the 4D harmonic 
oscillator Hamiltonian, is an easy task using eight dimension Clifford algebra. In addition, can be 
shown that using this Dirac operator on 4D-Moyal space, is possible define an eight-dimensional 
spectral triple. After defined the Dirac operator with the desired spectrum it is considered 
the total spectral triple as the tensor product of the "oscillating" spectral triple {A4,'H4,,T>4) 
with an almost-commutative triple and then is perform the previous described procedure of non- 
commutative geometry to compute the spectral action. We notice that matrix algebra introduces 
an extension of the standard potential in the commutative case, in fact the scalar field (j) and 
the Xfj^ fields are present together in a potential of the forrG0 (aX^ -k + /3(^* (p — 1)^, with 
a,l3 £ M+ and Xf^{x) = (0~-^)J^Xjy -|- A^{x) is a covariant coordinate. 

The high non-triviality of the vacuum makes very difficult to explicit the vacuum configura- 
tion of the system in [201 A. de Goursac, J.C. Wallet, and R. Wulkenhaar, using the matrix base 
formalism, have found an expressions from vacuum solutions deriving them from the relevant 
solutions of the equations of motion. Although, the complexity of the vacuum configuration 
makes the perturbative approach very complicated, in order to conduct some investigations will 
be considered a non-perturbative scheme using a discretized matrix model of the action in which 
the fields become matrices, the star product become the matrix multiplication and the integral 
turns in a matrix trace. 

Now comes in to play the numerical treatment, the standard method is to approximate the 
space by discrete points, for example using a lattice approximation and then calculate the ob- 
servables over that set of points [21]. Since an approximation in the position space is not suitable 
due to the oscillator factor of the Moyal product, instead the lattice approximation, will be used 
the matrix Moyal base, which was already used in the first renormalization proof of </?^-model 
restricted to finite matrices. Hence, will be performed a Monte Carlo simulation studying some 
statistical quantity such the energy density and specific heat varying the parameters Cl, ^^,Oi 
and gathering some informations on the various contributions of the fields to the action. The 
simulations are quite cumbersome due the complexity of the action and the number of indepen- 
dent matrices to handle but we are able to get an acceptable balance between the computation 
precision and the computation time. For the simulations is applied a standard Metropolis-Monte 
Carlo algorithm [22] with various estimators for the error and for the autocorrelation time of 
the samples. In general we chose the range of parameters in order to avoid problems with the 

^Einstein notation on repeated indices is used. 
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thermalization process, obtaining numerical simulations where is enough to wait a relative small 
number of Monte Carlo steps to compute independent results from the initial conditions, we are 
interested on the continuous limit that correspond to matrices of infinite size. We will consider 
various size of the matrices expecting a stabilization of the values of observables like the energy 
density, increasing the matrix size. In order to find same possible phase transitions will be used 
the specific heat which is a measure of the dispersion of the energy. The phase transitions are 
registered as peaks of the specific heat, increasing the matrices size. 

2 8-dim spectral action 

In this section will be computed a spectral action starting from a non-commutative spectral 
triple. The feature of this particular triple is the choice of a 4-dimension Harmonic Dirac 
operator. The idea behind this construction [18] is to relate the Dirac operator with the oscillator 
Hamiltonian operator. Roughly speaking, we look at the Dirac operator as a generalization of 
the Laplace operator so we have ~ H. Considering the spectrum a;(2n + 1) n G of 
the one-dimensional harmonic oscillator Hamiltonian H, can be deduced that is a non- 
commutative infinitesimal of order one. The non-commutative dimension of a spectral triple, 
equipped with the 4D Harmonic Dirac operator P| = i?4, is fixed by the non-commutative 
order of the inverse operator which is eight not four. This occurrence connects the spectral 
dimension to the phase space dimension instead the one of the configuration space [5] . In order 
to construct such Harmonic Dirac operator and the spectral triple we will work in the framework 
of the generalized n-dimensional harmonic operators. Will be studied the 4-dimensional case 
in order to construct the non-commutative spectral triple which is starting point for the field 
theory we are interested in. Having the 4-dimensional Harmonic Dirac operator with harmonic 
oscillator spectrum, to implement the Higgs mechanism we will consider the tensor product of 
the non-commutative triple with a finite Connes-Lott type spectral triple [11]]. We will fluctuate 
the total Dirac operator following the standard machinery [23^ [6] of non-commutative geometry 
to get "Gauged" Dirac operator. Thus we will proceed to compute the spectral action in which 
are present two U(l)-Moyal Yang-Mills fields unified with a complex Higgs field. 

3 Harmonic Dirac operators 



The Harmonic Dirac operator in n-dimensions can be defined using the Clifford algebra of 
M^" represented on the Hilbert space C^", it is very useful to consider n-dimensional fermionic 
annihilation and creation operators 6^, b\, and n-dimensional bosonic annihilation and creation 
operators a^, al satisfying for /x, = 1, • • • , n: 



Where = -^{ujx^ + d^), aji = -7={ujx^ — d^). Using this operators is possible to construct 
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a Dirac operator as: 

Vn = -iV2^6f"'al + iV2^6f"'a^ ^bl = i-^ {h^ + ftt ) + iuJx^' ® (6^ - b^) (3) 

summed over repeated index. We can define the fermionic part of the Hilbert space on which 
the Dirac operator ([3|) acts starting from the vacuum state by subsequent apphcations of the 
fermionic creation operators b\, on the vacuum 6|0) = 0, using the anti-commutation relations 
([2]) defining A(C"). The complete Hilbert space is Tin = S{W) A(C"). Beside, we can define 
a grading operator Xn as: 

n 

Xn = I'^llib^bl-b,bl) (4) 

Using the relations ([I])-© we can compute the square the Dirac operator ([3]) as: 

Vl = 2ujala^' O I - 2a;I O 6|,6^ = 2ujNb ®l-2ojlCdNF (5) 

Where Np and Nb are the number operators. In this form it easy to see that P^, being a "dif- 
ference" between fermionic and bosonic number operator, has only one zero mode corresponding 
to the vacuum state. For practical reasons it is convenient write as: 

Vl = u:6'"'{alau + a^aj.) I - 2a;I 5^"" {b^b^ - b^b^) = ® I + a; S„ (6) 

where in if„ we can recognize the harmonic oscillator Hamiltonian and the spin operator S„. The 
universality property of the Clifford algebra grants the existence of an isomorphism between the 
2n-dimensional Clifford algebra and the Hilbert space Tin = S {M."') C"^" . In this representation 
the Dirac operator is: 

P„ = ir'^9^ + a;r'^+% (7) 
Where F'^ turns to be F^ = (6^ + bji), F^+" = {{bj^ — b^) which satisfy the relations: 

FaFfe + FfeFa = 25ab with a, 6 = 1, • • • , 2n (8) 

Beside, the grading operator is represented as: 

Xn = {-ir{-l)^ ®Ti---T2n (9) 



4 An harmonic spectral triple for the Moyal plane 



In the framework of non-commutative field theories on 4-dimensional Moyal plane has been 
proved |24l [25] that the introduction of an harmonic oscillator term makes a i;^i^-model on 4- 
dimensional Moyal plane renormalizable. Such oscillator term can be written as: 



H„ 



+ Q'^Xi.x^ + 



(10) 



where x^ := 2(0 ^)^,yx'^, can be chosen as two copies of the Pauli matrix G = i9a2 I2 or 
explicitly: 

/ e o\ 

-6* 
61 
\ -9 J 



e 



(11) 
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With this choice we have ^ = — ^ (J5 12 . Quantum mechanics tell us that in the Hilbert space 
L^(IR^) exists an orthonormal basis ips, s G of eigenfunctions of Hm with eigenvalues 

A.M=(-(. + 2 + ^jj,.EN (12) 

The inverse extends to a selfadjoint compact operator on L^(R^) with eigenvalues \~^{rn). 
If we look at the trace the operator we find: 

00 

Tr(i/„^) = ^(n + 3)(n + 2)(n + l)(A„(m))^ (13) 

n=0 

Which is derived simply from the number of possibilities to express s as a sum of four ordered 
natural numbers. This means that belongs to the Dixmier trace ideal L(^'°°-'(L^(]R^)) of 
compact operators and the relation TP' = H implies that the 4-dimensional Moyal space has 
spectral dimension 8. From the previous section, we can define a proper Dirac operator just 
considering the 4-dimensional case obtaining a Dirac operator built from a 8-dimensional Cifford 
algebra: 

= iT^d^ + nv^+ix^, (14) 

Here, the Tk G Mi6(C), k = 1, 8 are the generators of the 8-dimensional real Clifford algebra, 
satisfying 

r^r; + TiVk = 25kii (15) 

We take the Hilbert space 7^4 = S''^(M^, S) = L^(M^) (g) C^^ of Schwartz functions of spinors over 
4-dimensional euclidean space. Accordingly with ([8]) for tjj G T-L^ we obtain: 

= ((-A + ^^x^.x^')! + S) V , S := -m{e-^Y^[T^, r^+4] (16) 

with A = d^d^. As algebra we chose the Moyal algebra Mg,: 

where (5(M^),*) is the algebra of the Schwartz functions on M^, with the Moyal product 

if^g){x) = J j d'y^J{x+-Q.k)g{x + y)e'^^^y) , f,geA, (18) 

The representation of the algebra on is by component-wise diagonal Moyal product [26] 
* : .44 X ^4 — )• 7^4. The Moyal product can be extended to constant functions using another 
representation of the product with the integral representation of the Dirac distribution. Taking 
in account, for smooth spinors, the identity 2x^^il: = x-ktp + ip-kx and the relation 

[x^,f], = iQ^Pdp (19) 

we compute the commutator of that action with the Dirac operator 



= iVt'iid^f) * + / * S^V') + \^T^'^\x^ ^ (/ ^ ^) + (/ ^ ^) ^ 3,^) 
-iWf * - ^^T^+Hf *(£;.* V) + / * (V- * x^)) 

= i(r^ + f]r^+4)(g^j)^^_ 



(20) 
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The previous commutator confirms that (^4, ^^4, P4) satisfy the mair0 axioms of spectral triple, 
in fact the commutator is bounded and due to its commutation with Moyal multiplication, order- 
one condition is fulfilled. Now we introduce a very useful relation connected to the heat kernel 
type expansion associated to a regular spectral triple taken from [19]. This relation will be used 
later in order to compute the spectral action. Considering a regular non-unital spectral triple 
{AjV-jV, J) and two pseudo-differential operator Aq £ ^'o(-^) ^1 £ of order respectively 

and 1. We consider the following decomposition: 

4 

^.tiV^+Ao+A,) ^ _ (21) 

j=0 

Using Duhamel principle [27] 

^-t(A+B) ^ ^-tA _i ^-st{A+B)j^^-{l-s)tA (22) 

Jo 

we can identify: 

E,it) = [ e-'''^\Ao + Ai)e-^''-''^'^' ■ ■ ■ {Ao + Ai)e 



and 



R = [ dsids2ds3ds4ds5e-'''^^'+^'>+^'\Ao + Ai)e-^'^-'''>'^' {Aq + Ai) 

J As 



(23) 



(24) 



The domains of the integrals Aj are the j-simplex: 

j 

Aj = {s G M^'; < si < S2 < si < • • • < < 1} ~ {s G M^+^; Si>0,^Si = 1} (25) 
Taking in account the relation: 

e-*^^^] = l^s-^^ (e-*^^'^e-*(i-^)^') = -t j\se-''^' [D^,A] e-'^'~'^^' (26) 

and considering the trace in [19] is computed the leading term of e-*(^^+^o+Ai) fQj- ^ — ). Q: 
Xr(e"*(^'+^o+^i)) 

= Tr({l - t{Ao + ^1) + ^(^0 + ^1^0 + AoAi + Aj) 

-j{Ao [V\ Ai] + Ai [V\ Ao] + Ai [V\ Ai] + AqA? + A^A^A^ + AIA^ 

+A?) + ^(^i [P2 [v\A,]]+2Al [V\A,]+A, [V\A,] A, + Ai)e''^'] 
+0{Vt) (27) 



Orientability axiom and Poincare duality will be not considered 
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5 4-dimensional harmonic Yang-Mills model 



Following the Connes-Lott models, in order to implement the Higgs mechanism, we consider the 
total spectral triple as the tensor product of the 8-dimensional spectral triple (^4, "^4, P4, Tg) 
with the two point Connes-Lott like spectral triple (C (E> C, C^, Mfii). The total Dirac operator 
of the product triple is: 

= P4 I + Tg O Mai (28) 



Or explicitly: 



Vr={ (29) 



The algebra becomes At = A^ © ^4 and acts by diagonal star multiplication ([TH]) on Ht = 
7^4 ® Hi. The fluctuated Dirac operator is found using Va = T>t + Siaj[^^T, bi] with Oj, ftj G At 
of the form {f,g), the computation of the commutator Dt with (/, 5) gives: 

L^{f)^ = f -kip is the left Moyal multiplication. From the commutator we deduce that the form 
of selfadjoint fluctuated Dirac has to be: 

^ ^/ v, + iT^ + nT^+,)LMn rgL,(<A) \ 

Where cp G Ai is the Higgs complex field and Af^,Bf^ £ A^ are real fields. The spectral action 
computation needs the square of '■ 



with 



D^(l) = dfj,<l)-iAf,-k(j) + i(j)-kBf, (33) 

Fa = {Vi, (P^ + ^T^+^)LMn} + (r^ + f^r^+4)(r. + nr,+4)L.(^^ * A'') 

= {LMn.id^ + ^^'M.(x^)} + (1 + f]2)L,(yl^ * A^') 

+ i Q[r^,r,] + io2[r^+4,r,+4] + f^r^r,+4) ^.(i^T)- (34) 

(M,(x^)?/^)(a;) = x^ijj{x) is ordinary pointwise multiplication and is obtained just replacing 

_ f) A _ fi A _ 

fill — '^fJ.-'^U '-'V-'^fl 



A with B. We can recognize in previous expression the field strength = duA^, — d^A 



^(^/^ * -^V ~ A/y * Afi) 

5.1 Spectral action 

Recalling the spectral action principle, the bosonic action can be defined exclusively by the 
spectrum of the Dirac operator. The general form for such bosonic action is: 

S{Va) = Tix{V\) (35) 
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Where x is a regularization function x '■ ~^ R+ foi' which trace exists. 
The trace in ([35]) is defined on B{L'^{R'^)) by 

Tr(^) = dx A{x,x) (36) 

together with the matrix trace including the Cliff'ord algebra. By Laplace transformation one 
has 

f'OO 

S{Va)= dtTr(e-*^i)x(i) (37) 



where x is the inverse Laplace transform of x{^)^ 

POO 

Xis) = / dte-^'m- (38) 
Jo 

The trace in ([37|) is given by: 

Tr(e-*^^)= / dxti{e-'^^{x,x) (39) 
Jr* 

Assuming the trace of the heat kernel e~ has an asymptotic expansion 

oo 

T¥(e-*^i) = ^ an{V\)r , <5 e N (40) 

n=—S 

we obtain replacing the previous expansion into ([3 



poo 

S{Va) = V an{V\) / em (41) 

To compute the integrals we have to consider separately the cases n ^ N and n G N: 

^ ^ / T^) ^-""-'Xis) forn^N 

\ (-1)"-V"H0) fornGN ^ ' 

Due to the nature of the x{t) function (usually one chose a characteristic function), we can 
assume x(0) much bigger then the derivatives X^^^O) for any m > appearing in (j42p . Conse- 
quently in the expansion (j40p we will take in account only the finite or singular part for i — t- 

Our strategy to compute the action is to use the relation ()27p . therefore after explicitly 
expressed and Ai we proceed to the calculus of the traces and in the end we will identify the 
leading part of the action comparing the result with the expansions (|40p - ()4ip . We can identify 
the operators and Ai appearing in the (j27p as follow: 

A.J-^^^-^')^:^^^^^^ _ ..V.,...J (44) 



2i{l + n^)L^{B>')V 



with 

VA,^ = (t)*^+{l + ^''){id^A^' + A^^A''), VB,^ = 4>^4> + {^ + ^''){id^B'' + B^^A^') (45) 
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V^^^ are define as v|f^^ = (p^ — iO^M,(5;^)) . We are allowed to split the traces in two 

parts a matrix trace and the continuous one. After the matrices trace computations we obtain 
[19] for the A field: 

Tr(e-*^^) 

= |l6cosh^(tn)tr(e-*-^4) - tr(16Vk,0) - t%{'i2i{l + Q^)A^) 

4.2 

+-T(iQVa,^ * Va,^ + 16(1 + n^)D^(t> * D^^^ + 8(1 + n^)F^,F^'' 

+2 , 

+-Tf, [32i{l + ^ + 32i(i + n^)VA,<p * A^' 

2 



-64(1 + n^fA" * a^^'^) + ^r^^ (-64(i + r?2)2^A' ^ 



t 



3 



g V ( - 64(1 + nyvA,^ * ^^A'^ - 64i(i + nyA" * a^i^^,^ 

+64(1 + n'^fAp * (5^^ A^^ + IdPdf'A") 

-64(1 + 1^2)2 (l^A,<^ ir Al" * A"" + Af" ir Va,^ ir A" + AI" * A" ir Va4 

-128i(l + {AP ★ (SpA'^) i^A^ + APir A^idpA") 



-jTpup (128(1 + n^ fAP * a''^^ - 128i(l + n^ fA" * * A'') 
+j%u (512i02(i + 02)(e~i)''M'^ * Ap) 

I- ^ I o2\4 yiu , ai'ap 4 0- , c-io-zi i r(2\3 



+^ Vpa ( - 256(1 + ^' fAP * &'dPA'' + 512i(l + ^' f A^ * ★ 



+256i(l + Q.^fA^ * (d^AP) * ^'^ + 256(1 + n^AP * A" * AP * A" 

+B field contribution + 0{Vi) (46) 
Where in order to simplify the notation we introduce the functions: 

T{f) = Tr,.2(M4) (L,(/)e-*^*) 
%,...p,{f) = TV^.(M4)(L.(/)V(^)---V£)e-*^^) 

V(/) = Tr^.(K4) (L.(/)v£)vWe-*^^) (47) 



The contributions for the B fields are obtained operating the following substitutions: 

{Ap ^ Bp, F^, ^ Fj^^, Va,4> ^ Vb,4>, Dp(t> o D^} (48) 

After the computation of the T,Tp^...p^, [§, we have all the ingredients required to compute the 
leading part of the action (j35p replacing all the traces into the (j27p . Using the trace property 
of the star product and the identities 

* / = \{~^,, f} + \ /] = \{^^. /} + (49) 
{Zf,J*g}^, = {zpJ}-kg-if-kdpg (50) 



'^Tiii' can be neglected 
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we get after some manipulations: 

TT{e-'^') = 2cosh\nt) 



where 



+ 2 V 2 6(1 + J]2)2 J (^A^- + * Fb ) 

+l±^{D^(t)*W4> + ^*D^'4>)\{z) + 0{Vt) (51) 



^,^ = T + ^M, ^.'^ = T + i?., Xl = i (52) 



Using the Laurent expansion of coth^(i') = t'^^ + + H + 0(t'^) and comparing the previous 
expression to the expansion (j4ip and putting xo = x(0) 'we are finally able to write the spectral 
action (1351) as: 



7r2(l + fi2)2 + ^2(1 + J^2)2 J ' 2 ' 6(1 + ^2)2 ) (^/^l * ^7 + ^/f. * 



2 



)2 ._x2 



"KtS^^°*^°"^) +2(1 + ^')^m'^^^^'a|(^) + 0(xi) (53) 

We notice that Higgs mechanism introduces an extension of the standard Higgs potential in the 
commutative case, in fact the Higgs scalar field (p and the X'^, fields are present together in 
the potential. In this way the gauge field takes part in the definition of the vacuum. Another 
important property of the action, considering the X^, X^ as independent, is the invariance 
under the translations: 

0(x) ^ 0(x + a), X'Xix) XfXix + a), X^(x) ^ X;^(x + a), X^{x) X^{x + a) (54) 

which in other (/>^-renormalizable theory is broken. Beside, the action is invariant under U{1) x 
[/(I) transformations: 

(p ^ UA-*: 4>-*:Ub, X ^ ua* X'^-kUA, X'^ ^ ub * X'^ -kUB (55) 
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In field theory the ground state can be defined through the minimum of the action, the relevant 
part of the ()53p for the minimization is: 

-2 (lf^*o * " + + a-')D,4,*D^\(z) + 0(xi) (56) 

Where we have omitted the constant part and we have rescaled the coefficient in front of the 
integral. Considering the fields X^, as fields variables instead A^^, we can state that 
each terms of the action is semi-positive defined, so in order to find the minimum it is sufficient 
to minimize them separately. There are the two possible minimum for the field strength part 
and for the covariant derivative part: the trivial solution with cp and X^, equal to the null 
fields and the solution with (j), X^, proportional to the identity. In each cases both the field 
strength part and the covariant derivative part disappear. For the potential parts we have: 

for (/> = X^ = X^ = 
for = a^I, X^ = uaI^, X^ = asP 



Referring to the second case and minimizing the potentials, the minimum seems to be for 

<^o = ^ cos al, X% = X^o = ^ ^J^^JI^l sin al^ (57) 

However, the previous position is not allowed because the identity does not belong to the algebra 
under consideration. In general the non-triviality of the vacuum makes very difficult to explicit 
the vacuum configuration of the system in [20] A. de Goursac, J.C. Wallet, and R. Wulkenhaar, 
using the matrix base formalism, have found an expressions from vacuum solutions deriving 
them from the relevant solutions equations of motion. Although, the complexity of the vacuum 
configuration makes the perturbative approach very complicated, in order to conduct some 
investigation in the next section will be consider a non-perturbative approach using a discretized 
matrix model of the action ()56p obtained using a Moyal base. In this setting the action reduces 



Va = Vb 



VA = (ai + 4 



2 , . 2 



X-i 
Xo 

Xo 
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to 



S{Va) = 

I 1 / .4 J ( ji-^r _ (1 - \ ^ , ^ j,,u. 

+ (1 + 1^2)2 y ^ 1^ 2 6(1 + 1)2)2^ + i 

+2(l + f]2)i)^^^;D^|(^)+0(xi) (58) 

The omitted factor for the finite matrix model of size N becomes constant so can be ignored. 
The minimum is obtained hke before and formally is ()57p in this case the identity, of course, 
belongs to the matrix space. It is interesting to notice that the vacuum of the finite model, due 
to the Higgs field, is no longer invariant under the transformations ()55p . but is invariant under 
a subgroup of U{N) x U{N): 

UA = ub — > UA*i* ub = ua* ua = I (59) 

Having discretized the model will be performed a Monte Carlo simulation studying some statis- 
tical quantity such the energy density, specific heat, varying the parameters Q, o and gath- 
ering some informations on the various contributions of the fields to the action. The simulations 
are quite cumbersome due to the complexity of the action and to the number of independent 
matrix to handle. 



6 Discretization of the action 

The first step across the numerical analysis is to apply a discretization scheme. Various schemes 
can be used like lattice approximation, but the nature of star product due to its oscillator 
exponential, makes the lattice approach not suitable without adaptations. We will use another 
discretization scheme in which our fields are approximated by finite matrices and the star product 
becomes the standard matrix multiplication. Using the identity Dfj_(j) = (p-k Xb^ — Xa^ * </> we 
can recast the action (|58p in the following form: 

+2(1 + n^) * Xs^ - * * 1^ - * I (z) 

+Oixi) (60) 

As a first approach to the numerical simulation and forced by limited computation resource, 
we will consider the Monte Carlo simulation of the previous action around its minimum and 
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the simulation will take y as a positive parameter. In this setting the behavior of the 
simulations will be identical for the negative case and avoiding to be negative we have 

/CO 

not any problems about the thermalization. In order to define the previous action around the 
minimum we translate the fields </>, X^^, XB^l using the following translated fields: 



il) + . j^—^ cos al (61) 
Xo 



Xb, = + (63) 



Substituting the previous fields into (j60p we get a positive action with minimum in zero: 

+ [il)-ki) + ji COS a{il) + + CY^ ★ Ya^, + ^ll^YA^, sin a) ^ 
+ ill) -k il) + ^l cos a{'il) + ^) + CYj^ ★ Yb^ + /iF Yb^ sin a) ^ 

+2(1 + n^) {(Yb^, - YA^,)fi cos a + * Yb^ - Ya^, * V) 

*((y-^-y^)^cosa + v^*yj'-y^*V^) |(z) + o(xi) (64) 

Where for simplicity we put: 



1+"" „ _ (i-gf (i-n'j* x-i 



6.1 Discretization by Moyal base 

The following treatment is mainly taken from j26^ [28] as introduction to the Moyal base which 
will be used later. We can define on the algebra ffi|, a natural basis of eigenfunctions fmn of the 
harmonic oscillator, where m, n G N. This base satisfy the ★-multiplication rule: 

ifmn*fkl){x) = 6nkfml{x) (66) 

and this useful property: 

d^xfmn{x) = Smn I d^xfo = 27r6'(5mn (67) 



The previous multiplication rule associates the ★-product between f^i with the ordinary matrix 

^1 



product: In this base we can write any elements of Rq but we have to require the rapid decay 



[28] of the sequences of coefficients {amn}- 

OO CO 1 

<^mnfmn{x) G M|, if, ^ ((2m + l)2^(2n + lf^\a^n?) ' < OO for ah k. (68) 



m,n=0 m,n=0 
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The eigenfunctions fnm can be expressed with the help of Laguerre functions [251 IMl EI] '■ 



n—m 



Our fields can be expanded in this base as: 

X^{x) = V] X'i^n-t^ frnini{xo,Xi)fm2n2{x2,X3) (70) 

7712 ^2 

and 

^{x) = V'™i"i/mini(2;o,a;i)/m2n2(2;2,2;3) (71) 

Using this base we can forget the Moyal product in this way the model becomes to 9-matrix 
model. A ^-product between two fields using ([66]) can be written as 

^'(x)*^>(x) = ^' mini $fci!i/mini (3:0, Xi)*/fe,ii(xo,Xi) 

?Tii,ni,fci,iiSN 
X fm2n2ix2,X3) -k fk2l2ix2,X3) 

= y ^"'i"i^kj^h^nikiSn2k2fmih{xo,Xi)fm2l2{x2,X3) 

^ — ' "^2"2 k2l2 

mi,ni, fci,iiGN 

= 'I'^>miii/„,/,(X0,Xl)/„2«2 (2^2,3:3) (72) 



where 



*$mi;i = *mini$„^,^ (73) 



m2;2 ™2"2 n,; 

™i,n-2£N 



So the star product became a "double" matrix multiplication, the action, the equations of field 
and all treatments can be conducted on the infinite matrices instead directly on the continues 
fields. 

Beside, for finite matrices, the N^-indexed double sequences can be written as tensor products 
of ordinary matrices, 

K 

X™ini='S^X* (74) 

1=1 

Since the matrix product and trace also factor into these independent components, the action fac- 
tors into S = S{iJ^\Yl\ Y^')S{'il^'^\Yl\Y^'). Then, regarding ah i;^\Yl\ Y^, i?\Yf,Yl' 
as random variables over which to integrate in the partition function, the partition function fac- 
tors, too: 

|p(V'^\yi\yi\V'2\yr,yJ^)...p(^i^,yi^,yr,^2^,yr,yr) e"^ 

= (I p(V'^\yi\yi\ v'^yi^yf ) e-^'^^'^y): y'B)-s{^'\yJ: yl')^ . (75) 

We may therefore restrict ourselves to -ftT = 1. Using this approximations the calculus will 
be performed just on standard infinite matrix, but is not enough to be handled numerically. 
We have to perform a truncation in order to obtain finite matrices, this truncation will consist 



15 



in a maximum m,n < N in the expansion (|71|) -(|7U|). It is easy to verify that this kind of 
approximation corresponds in a cut in energy, in fact from the definition of fmn we have: 

{H,fjnn}* = 0{m + n + ^)fmn{f,P) (76) 

Beside, can be proved [31] that the functions fmn with m,n < N induce a cut-off in position 
space and momentum space: 

Pmax ~ V29 , for m, n < (77) 

and 



Pmax ^ y-j- for m, n < iV (78) 
Summarizing, to operate the discretization we have the following correspondences: 

0(x) G Rq ^ G Mjv (79) 

Y^{x) G R| ^ G Mtv (80) 

Y,^ix)eR% ^ y^^(^M^ (81) 

a{x)dx Tr(a) (82) 

After truncating the representative matrices is convenient to operate another substitution |20] : 





— -^0 




Zo 


- Y^ 


-iY^ 


Zi 


- Y^ 


+ iYf, 


Zi 


- Y^ 


-iY^ 


Z2 


- Y^ 


+ iY2^, 


Z2 


- Y^ 




Z3 


- Y^ 

— ^2 




Zz 


- Y^ 

— ■'2 





(83) 



In the end the discretized action is: 



With 



£'4 = ^ _^ ^2 '^^ + "^^0 + ^Y + ^DqC^Do + Cd^Cdx + Cd2^D2 + ^DsJ-^Ds) (84) 



CiF = ^([Zo,Zo]'+ [^l,^l]' + i([^0+^0,^2-^2]'- [^0 + ^0,^2 + ^2] 

+ [Zo~Zo,Z2 + Z2f - [Zo-Zo,Z2~Z2]^ - [Zi+Z^,Z3 + Z3]^ 

+ [Zi + Zi,Z3~Z3f + [Zi-Zi,Z-s + Z3]^ - [Zi-Zi,Z3-Z3]^^) 

C4V0 = (i/"i/; + MCOS 0(1/- + 1/.) + i ({^0,^0} + {^2,^2}) 

+ ^^((-1 + 0(^0 + Z2) + (1 + ^)(Zo + Z2)))' 

£4Vi = {^PtP + ^lcosa{tP + ^P) + ^{{Zl,Zl} + {Z3,Z3}) 

+ ^^((-1 + 0(^1 + Zs) + (1 + ^){Z^ + Z3)))' 

CiDo = \/2(l + r!2) cos q(^i +Zi- Zq - Zo) + V(Zl + Zi) - (Zo + Zo)t/') 

£4Di = \/2(l + n2) COS q(Zi - Zi - Zq + Zo) + V(Zi - Zi) - (Zo - Z„)V') 

£4D2 = \/2(l+n2) (^ COS a(Z3 + Z3 - Z2 - Z2) + V(Z3 + Z3) - (Z2 + Z2)V') 

£403 = ^/2(1+W) {ficosaiZs -Z1-Z2+ Z2) + tpiZa - Z3) - (Z2 - Z2)V') 
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Where for simplicity we have omitted the hat on the matrices and the bars stand for the hermitian 
conjugate. In this case the action (j84p becomes 5 complex matrix model instead eight real 
matrices Ya^, Yb^ and one complex matrix ip. This form may seem cumbersome but it is more 
comfortable for numerical simulations. The next step is to define the estimator for the average 
values of interest and to develop some numerical parameters in order to analyze the numerical 
results. 



7 Definition of the observables 

Following Monte Carlo methods, will be produced a sequence of configurations {{tp, Zi)j},j = 
1, 2, • • • , Tmc and evaluated the average of the observables over that set of configurations. The 
sequences of configurations obtained, a Monte Carlo chain, are representations of the config- 
uration space at the given parameters. In this frame the expectation value is approximated 
as 

(0)=»^i:0, (85) 

where Oj is the value of the observable O evaluated in the j-sampled configuration, {■ilj,Zi)j, 
Oi = OKip, Zi)j]. The internal energy is defined as: 

E{Q,fi,a) = {S) (86) 

and the specific heat takes the form 

C{n,f,,a) = {S'')-{Sf (87) 

These terms correspond to the usual definitions for energy 

1 dZ 

and specific heat 

dE 

C{n,fi,a) = — 

where Z is the partition function. It is very useful to compute separately the average values of 
the four contributions: 

F{n,fx,a) = {Tr Cf) (90) 
Voin, fi,a) = (Tr£vb> (91) 
X^i(0,/i,a) = {TvCv,) (92) 
D{n,^i,a) = {Tv{CD,Cr).)) (93) 
V = Vo + Vi (94) 

Where i,j = 1, 2 or i, j = 1, • • • ,4. 
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7.1 Order parameters 



The previous quantities are not enough if we want to measure the various contributions of dif- 
ferent modes of the fields to the configuration ^/^, Zj. Therefore we need some control parameters 
usually called order parameters. As a first idea we can think about a quantity related to the 
norms of the fields for example the sums Ylnm iV'nmP, X^nm l-^mmP of all squared entries of our 
matrices, this quantity is called the full power of the field [291 ED] and it can be computed as 
the trace of the square: 

= Tr(|^|2) (95) 
Zl = Tr(|Z,|2) (96) 

However, {(pa) alone is not a real order parameter because does not distinguish contributions 
from the different modes but we can use it as a reference to define the quantities: 

Af 

2 _ I |2 

^0 — / - \(^nn 



E 

n=0 

TV 



n=0 

Referring to the base ([69]) it is easy to see that such parameters are connected to the pure 
spherical contribution. This quantity will be used to analyze the spherical contribution to the 
full power of the field. We can generalize the previous quantity defining some parameters (pi in 
such a way they form a decomposition of the full power of the fields. 

^l = vl + Y.'^l Zl = Zl + Y,Zl (98) 
I I 

Following this prescription the other quantity for I > can be defined as: 

I I 

= \anm{l - 6nm)\'^ , Zfi = ^ | Zjnm (1 " 5„m) | ^ (99) 

n,m=0 n,m=0 

If the contribution is dominated from the spherical symmetric parameter we expect to have 
i^a) ~ (<^o)' (Zia) ~ (Zio)- With (j99]) we can define the order parameters (fl, Zf^ as a particular 
case / = 1 we have 

"Pi = kioP + |aoiP, Z\ = \ziXQ\' + l^ioil^ (100) 

In the next simulations we be evaluated the quantities related to / = and to Z = 1 as represen- 
tative of those contribution where the rotational symmetry is broken. Using higher I in (j99p we 
can analyze the contributions of the remaining modes, turns out that the measurements of the 
first two modes are enough to characterize the behavior of the system. 



8 Numerical results 

Now we discuss the results of the Monte Carlo simulation on the approximated spectral model. 
As a first approach we use some restrictions on the parameters. Will be considered the spectral 
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action around its minimum (j84|) in order to simplify the calculation, in this frame in the action 
is symmetric under the transformation ^ — )• —ji so will be used the condition /i > and /x^ > 0. 
Will be explored the range fi G [0, 3.1], this interval was chosen to show a particular behavior of 
the system for fixed il.. The parameter il. appears only with its square and is defined as a real 
parameter, therefore for the $7 too we require > 0. Beside, if we refer to the scalar model, 
is possible to prove [23] using the Langmann-Szabo duality [32J, that the model can be fully 
described varying in the range [0, 1], for higher the system can be remap inside the previous 
interval. In the present model the L-S duality does not hold any more, but forced by limited 
resource, we conjecture that the interval Q € [0, 1] still enough to describe the system. The last 
parameter to consider is a, it is connected to the choice of the vacuum state, the range of a is 
[0,27r]. The study of the system varying this parameter is quite important from a theoretical 
point of view because is related to the vacuum invariance. Beside, in the action appear some 
contributions proportional to {sma)/Q which seem to diverge for = 0, numerically we have 
verified that is an eliminable divergence and the curves of the observables can be extended in 
zero by continuity. Studying the dependence on a we can conclude that in the limit N ^ oo 
the observables are independent from a, therefore for our purposes a will be fixed equal to zero 
avoiding the annoying terms. In general for each observable are computed the graphs for N =5, 
10, 15, 20 matrix size. 

8.1 Varying a 

We start looking at the variation of energy density and the full power of the fields density for 
fixed /i and $7, varying a £ [0, 2tt]. As representative here will be presented the graphs for ^ = 1, 
= {1, 0.5} but we obtain the same behaviors for any other choice of the parameters allowed in 
the considered range. All tree graphs show an oscillating behavior of the values, this oscillation 




Figure 1: Total energy density and full power of the fields density for ip^, Zq^ (from the left to the right) fixing 
fi = 1, = 1, varying a and N. N — 5 (circle), A'^ = 10 (square), A'' = 15 (triangle), A'^ = 20 (cross). 

is present in all other quantities measured. The amplitude of this oscillation becomes smaller 
and smaller increasing the size of the matrix and this is true for all the quantities measured. 
The same trend is described in figl2] in which position of the maximum are different but the 
amplitudes becomes smaller increasing This results allow us to consider a = for all next 
graphs, since we are interested in the behavior of the system for A — t- oo. This occurrence 
simplify all the next simulations thanks to the vanishing the of terms ~ (sin a)/^} appearing in 
the discretized action. Beside, such results induce us to reckon the parameter a as connected to 
the remaining invariance of the vacuum state for the exact model. 



19 



Figure 2: Total energy density and full power of the fields density for (^^, Z^^ (from the left to the right) fixing 
/I = 1, n = 0.5, varying a and N . 

8.2 Varying Q 

Now we will analyze three cases in which /i is fixed to 0, 1,3, a is zero and we vary S [0, 1]. 
In the rest of this section we ignore for the computation of D, F the prefactor (1 + fi^). In 
this way we focus our attention to the integral as the source of possible phase transitions. The 
graphs in figj3]show the total energy density (5)/7V2 ^nd the various contributions: the potential 
V/N'^, the Yang-Mills part F/N'^ and the covariant derivative part D/N"^, for = 1. There 
is no evident discontinuity or peak and increasing the size of the matrices the curves remain 
smooth. Comparing the energy density and the various contributions figS] we notice that the 




Figure 3: Total energy density and the various contributions for fi — 1, a = varying Q and A'^. From the left 
to the right E, V, D, F with iV = 5 (circle), = 10 (square), iV = 15 (triangle), iV = 20 (cross). 

contributions between F and V balance each other and the total energy follows the slope of D, 
this behavior continues increasing the size of the matrices figSl 
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Figure 4: Comparison of the total energy density and the various contributions for ^ = 1, a = 0. E (circle), F 
(triangle), D (cross), V (square). With iV = 5 (left) and iV = 20 (right). 




Figure 5: Specific heat for ^ = 1. 

The specific heat density shows figl5] a peak in = 1, this peak increases as N increase. 
This behavior is typical of a phase transition, the peak is not clear for small A'^ due to the finite 
volume effect. In order to gain some informations on the composition of the fields we look at 
the order parameters defined in the previous chapter. Starting from ip field, in the figure [5] it is 
showed the graphs for 99^, (/?q and (p\ for N = 5. 

The three values 99^, (/?q and ipf seem essentially constant, comparing the three graphs figE]it 
is easy to see the dominance of the spherical contribution cpQ to the full power of the field. The 
behavior of the Zq fields is different, the spherical contribution becomes dominant approaching 
to = starting from a zone in which the contribution of Zqq and Zq^ are comparable. For 
brevity will not be showed only the graphs for Zq^, Zqq and Z^^ but taking in account the 
statistical errors the other Zi related graphs appear compatible to the Zq case. The dependence 
of the previous quantities on N are showed in the following graphs figlZl The values of the 
quantity of all the previous parameters decreases with but the dominance of the on the 
total power of the field is independent by N. The peak related to Zq decrease, but if look at 
the single graph for N = 20 the spherical contribution approaching the point J7 = it features 
a peak. 
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Figure 6: On the left comparison of ip^ (circle), ip^ (square) and ipl (triangle) density. On the right comparison 
of (circle), Zqo (square) and Zqi (triangle) density. 




Figure 7: Starting from the up left corner and from the left to the right the densities for </3o, Vii ^Oa, Zqq 
and Zqi for fJ. — I varying Q, and A^. 

Now we will analyze the model for n = 0; figE] shows the graphs for total energy density and 
the contributions V, D, F. The slope of the total energy density seems to be constant. The 
D contribution and the F do not balance each other like in the previous case, but all the three 
contributions balance among them self to produce a constant sum. 
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Figure 8: Total energy density, various contributions and the comparison among them for /i = varying f2 and 
A*'. From the left to the right E, V, D, F and comparison. 





Fi gure 9: Specific heat density for fi = varying n and TV. 



The specific heat density shows figE] again the peak in J7 = 1 as A'^ increase. For the other 
quantities ip"^, ip^, (p\ and Zq^, Zqq, Zq^ we have the same behavior fig(12]of = 1 case, except 
for the oscihation appearing in the Zq^, Zqq graphs close to zero, anyway it appears only for 

= 5. A complete different response of the system is described in the graphs for /i = 3, as 
we can see from figllll The slope of total energy density is very similar to the F component 
instead D. Beside, appears a maximum around = 0.4 for A^ — )• oo. This dramatic change of 
the graphs might be interpreted as consequence of a phase transition ruled by the parameter ^, 
actually in the next section we will find a peak in the specific heat density for some fixed Vt and 
varying ^ S [0, 3]. 
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Figure 10: Starting from the up left corner and from the left to the right the densities for ip\, cpg, (pi, Zg^, Zqq 
and Zqi for fi = varying Q and N. 




02 U4 tli> OS la 



Figure 11: Total energy density, various contributions and the comparison among them for ^ = 3 varying D, 
and A'^. From the left to the right E, V, D, F and comparison. 

Specific heat density displays figll2]a strong change too, in fact instead the peak in = 1, 
it appears in the opposite side of the studied interval in J] = 0. This peak too, due to its grows 
increasing N could indicate a phase transition. The figll3l describes the behavior of the order 
parameters densities ip^, ^Pq, ^pf and Zq^, Zqq, Zq^, they have a similar aspect to the previous 
relative graphs. For the V field the spherical contribution remains dominant, beside in the iff 
graph appears a deviation from the constant slope this deviation is evident for N = 5 but still 
there for higher A^. The order parameters for Zq display a peak close to the origin without 
oscillations even for N = 5. This maximum for higher does not move closer to the origin, 
in other words, this shift in not due to the finite volume effect. Even for Zq^ graph appears a 
deviation from the constant slope, a small peak which becomes shifted and smoother for higher 
N 
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Figure 12: Specific heat density for jj, — 3 varying fl and A'^. 




Figure 13: Starting from the up left corner and from the left to the right the densities for (^9^, c/Pqi V-'i? ^Oa^ ^oo 
and Zqi for = 3 varying Q and A''. 

8.3 Varying 

In this section is analyzed the response of the system varying /i G [0, 3] where is fixed to 
0,0.5, 1 and a is always zero. We start displaying the graphs figJ14l of the total energy density 
and of various contributions for = 0. There is no evident discontinuity but appears a peak in 
the total energy density around /i ~ 2.5 for = 20. Comparing all the contributions is easy to 
notice that the slope of the total energy is dictated by the curve V of the potential part. 

As mentioned before, the specific heat density figJTS] features a peak around /i ~ 2.5 for 
= 20 and again, due to this behavior as N increase, we could relate this peak to a phase 
transition. The plots for the quantities ip'^ and denote a strong dependence on /i, in particular 
the slope of (^q seems mostly linear, iff related graphs feature a similar behavior but the slope 
is no longer linear. Comparing the three graphs figdU] we deduce that close to the origin the 
non spherical contribution is bigger the spherical one, increasing fi this situation capsizes and 
(/?Q becomes dominant respect ipf. The behavior of the Zq fields is quite different, referring to 
figure \T6\ the spherical contribution is always dominant for the all interval fi £ [0, 3]. The curves 
for Zq^, Zqq are compatible to the constant slope, for Z^^ we have the same dependence on /i 
in particular there is a smooth descending step, however this step becomes smoother for bigger 
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Figure 14: The total energy density and the various contributions for SI = varying /i and A^. From the left to 
the right E, V, D, F an comparison with N = 5 (circle), TV = 10 (square), A*' = 15 (triangle), A'^ = 20 (cross). 
For the comparison: E (circle), V (square), D (triangle), F (cross). 



F 
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Figure 15: Specific heat density for Q = varying /i and A'^. 

N. It is behooves to say that due to some cancehations effects the statistical errors are quite 
big and they can hide some dependence, anyway this results tell us about the dependence of the 
order parameter for Zi and in general of the system, on the two choice = or 7^ 0. 

Now we will analyze the model for 0, = 0.5; as figJ 171 shows the graphs have a different slope 
comparing to the previous case, the maximum of total energy density follow the one of the V 
component. If we focus ourself only on the total energy graph and we compare it with the one 
for = 0, we notice a shift of the maximum for each A^. In particular in figllTI some maximum 
are moved outside the considered interval. We can find this shift very clearly looking at specific 
heat density graph figUHl we find again the peak as increase but it is shifted around fi ^ 3.3. 
The graphs figJ21lfor </j^, ip^ have the same behavior of O = 0.5 case, excluding some fluctuations 
close to the origin for N = 5 due to the finite volume effect, ipl graph displays an almost constant 
curve. However, close to the origin, the spherical contribution and the first non spherical one 
are comparable. The introduction of 7^ creates, in the Zq fields order parameters fig{T9l a 
dependence similar to the graphs for the ip; the full power of the field density and the spherical 
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Figure 16: Starting from the up left corner and from the left to the right the densities for i^^, Voi Vii ^Oa, ^oo 
and Zqi for f2 = varying fi and TV. 




Figure 17: Total energy density and contributions for D, = 0.5 varying /i and A'^. From the left to the right E, 
V, D, F. 

contribution are no more constant and they grow increasing fi. Even in this case the spherical 
contribution is always dominant excluding the region around = 0. 
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Fi gurc 18: Specific heat density for ^1 = 0.5 varying fi and TV. 




Figure 19: Starting from the up left corner and from the left to the right the densities for (p^, ^po, fi, -^oa, Zq: 
and Zqi for SI = 0.5 varying p and A'^. 



The last set of graphs for the 4-dimensional model are obtained fixing i7 = 1, due to the 
vanishing of prefactor in front of the Yang-Mills part of the action the F contribution is always 
zero. The following diagrams for the energy and contributions show the absence of the previous 
peak and comparing again them with former graphs they seem a sort dilatation. The specific 
heat density does not show the peak in zero any more figl2I] and the curves does not show any 
particular point as increase, actually the peak can be found for higher ^. At last in figi22 
we found a behavior of the order parameters density for Zq and -0 fields similar to the former 
graphs for Vt = 0.5 and they are compatible with a dilatation of the previous diagrams. 
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Fi gur6 20; Total energy density and the various contributions for = 1 varying /i and A^. From the left to the 
right E, V, D. 
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Fi gurc 21: Specific heat density for = 1 varying /i and A^. 

Conclusions and prospectives 

We have presented a study of a spectral action model constructed in order to extend to the Yang- 
Mills theories the non-commutative Wulkenhaar-Grosse model. The main aim of this work was 
to test a first Monte Carlo approach based on a non-perturbative regularization method. We 
have performed Monte Carlo simulations and obtained the values of the defined observables 
varying the parameters of the system. Despite the complexity of the approximated spectral 
action considered here we were able to obtain some reliable numerical results, we can conclude 
that a numerical approach of this kind of model using the matrix Moyal base seems feasible. The 
specific heat density shows various peaks indicating phase transitions, in particular studying the 
behaviors for some fixed /j, we found a peak around = for fi = {0, 1} and peak in = 1 for 
/i = 3, beside we notice a huge change in the energy density and in its contributions between 
the cases = {0, 1} and /U = 3. Other peaks in the specific heat density can be found varying 
fi and fixing 0, the graphs show that increasing the peak in the specific heat, starting from 
/i « 2.4 for fi = 0, is moved towards higher fi. The order parameters introduced show a strong 
dependence on the occurrence of fi = or fi 7^ 0. Referring to the fixed fj, graphs we found 
a peak in the spherical contribution for the gauge fields Zi, we can interpreted this slope as a 
sort of symmetry breaking introduced by fi 7^ 0. Additionally, varying fi and fixing fi the other 
parameters display an increasing slope with fj, for all fields and all situation but one; the graphs 
of the order parameters concerning Zoa, Zqq for fi = show a constant behavior. The natural 
next steps in the numerical study of this model, could be the computation of the transition 
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Figure 22: Starting from the up left corner and from the left to the right the densities for ip\, cpg, (fii, Zg^, Zqq 
and Zqi for f2 = 1 varying n and N. 

curves in order to separate the phase regions and classify them using eventually some additional 
order parameters. Our treatment, forced by limited resource, was conducted conjecturing that 
the system can fully described varying Q in the range [0, 1] but since the L-S duality does not 
hold any more in our case, will be very useful to extend this range. Actually, the computed 
graphs does not show any periodicity in 17 G [0, 1] so there are some hints to infer that the 
range [0, 1] is not enough, however only a direct computation will clarify this point. Will be 
very interesting to do not require any more the condition /i^ > 0, in order to implement this 
change we have to conduct the calculation no more around the minimum of the actions. The 
expansion of the parameters space, together with the classification of the different phase regions 
allows us to compare our model with the results of the simulation conducted on the fuzzy spaces, 
looking in particular about the occurrence of the so called non-uniformly ordered phase which 
is connected with the UV/IR mixing. Since we have constructed our model starting from a 
renormalizable one this study is very desirable. 
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